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Abstract 

The motion of the structure determining components is highly collective, 
both in amorphous solids and in undercooled liquids. This has been deduced 
from experimental low temperature data in the tunneling regime as well as 
from the vanishing isotope effect in diffusion in glasses and undercooled liq- 
uids. In molecular dynamics simulations of glasses one observes that both low 
frequency resonant vibrations and atomic jumps are centered on more than 
10 atoms which, in densely packed materials, form chainlike structures. With 
increasing temperature the number of atoms jumping collectively increases. 
These chains of collectively jumping atoms are also seen in undercooled liq- 
uids. Collectivity only vanishes at higher temperatures. This collectivity is in- 
timately related to the dynamic heterogeneity which causes a non-Gaussianity 
of the atomic displacements. 

I. INTRODUCTION 

Research into the physics of glasses spans a large temperature range, from the liquid to 
near K. There are a number of characteristic temperatures. The melting temperature, 
T m , marks the crossover to the undercooled regime, where the liquid is no longer the ther- 
modynamic ground state. The Kauzmann temperature, Tk, is the temperature where the 
entropy of the liquid would fall below the one of the corresponding crystal. It is sometimes 
interpreted as an ideal glass transition temperature. Before Tk is reached upon cooling, any 
real system arrests at the glass-transition temperature, T g , into a solid like state, the glass. 
This temperature T g shows e.g. as a change in volume expansion or a (smeared) jump in 
the specific heat. By definition one has in all systems Tk < T g < T m . The exact value of the 
glass-transition temperature depends on cooling rate, pressure and other parameters Q. 

In many materials one observes a sharp upturn of the viscosity curves at temperatures 
well above T g . This change in the dynamics of the melt is taken in theories such as the mode 
coupling theory (MCT) || as the true transition which defines a temperature T c > T g . 

At the lowest temperatures, below w 1 K, the dynamics in glasses is dominated by two- 
level systems [|[ which can be described by the tunneling model which was formulated nearly 
thirty years ago In addition to these tunneling states, one observes local relaxations 

in glasses and at frequencies near 1 THz an excess of vibrations compared to the Debye 
spectrum, given by the sound waves. This excess leads to the "boson peak" in the inelastic 
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scattering intensity. To describe the excess low energy excitations the tunneling model was 
extended to the soft potential model |],f7| . From this model "universal" temperature depen- 
dencies can be derived for temperatures of some 10 K. Comparing the model to experiment 
one finds that the excess low energy excitations are collective motions of 10 to 100 atoms 
H|9|]. Similar values were found by computer simulations for the low frequency resonant 
modes |T0|| . 

The most popular model to describe the under-cooled liquid is the mode coupling theory 
0. In this statistical model a nonlinear memory kernel leads to a blocking of modes at T c . 
It predicts scaling relations and dependencies of the form (T — T c ) 7 for quantities such as 
the diffusion constant. Many aspects of this theory have been at least semi quantitatively 
verified by computer simulations, see e.g. [[Tl|1 . 

At present there is no common theory for both the glassy and the super-cooled state. 
The transition from one regime into the other is probed by experiment and increasingly 
by computer simulations. The accumulated data should help to probe the limits of the 
existing theories and guide towards a theory bridging the present gap between glass and 
liquid theories. 

In this contribution we will review computer simulations both in the glass and in the 
liquid. The emphasis will be on collective effects observed in both. 



II. SIMULATION MODELS 

Our molecular dynamics simulations are done for systems of 500 to 10000 atoms with 
periodic boundary conditions. The equations of motion are integrated by the velocity Verlet 
algorithm with a time-step of order fs. We use cooling rates of > 10 11 K/s and aging times 
of order ns. For the model systems of soft spheres and Lennard- Jones atoms these values 
correspond to a conversion to Ar. Zero external pressure was exerted on all but the soft 
sphere systems where the volume was kept constant. 

To check the material dependence of the results we use model systems (soft spheres (SSG) 
and Lennard- Jones (LJ), both monatomic and binary) and models of the binary metallic 
glass CuZr and of Se. Details of the simulations can be found in the references given further 
down. 

The soft spheres (SSG) are described by a purely repulsive potential 

V(R)=e(a/R) n + V cut (1) 

where n = 6 or n = 12 and V cn t is a small correction to give a smooth cutoff. 

As a simple model potential which allows for zero external pressure we use a Lennard- 
Jones (LJ) potential 

V(R) = 46 [{a/Rf 2 - (cr/i?) 6 ] + V cnt . (2) 

We study both monatomic and binary LJ-systems. In the latter case we use the parameters 
of Kob and Anderson ]TT|] with a different cutoff. 



As a typical example of a binary metallic glass we simulate Cu33Zr 67 , described by a 



modified embedded atom potential. For details see [12]. The results do not strongly depend 
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on the choice of interaction. As far as data are available there is good agreement with the 
similar NiZr system modeled with a totally different interaction 

As an example of a good monatomic glass former we choose selenium. Se has a coordina- 
tion number near 2. The amorphous structure consists of inter-netted rings and chains. We 
describe it by a potential of the Stillinger- Weber type. In this potential the covalent bonds 
are described by an additional three body term V3 (|R X — R 2 |, |R X — R 3 |, cos 6213) where 
O213 is the apex angle [I3|. The parameters were fitted to both molecular and crystalline 



data. 

The examples presented in the following will be taken from simulations of these different 
materials. We will emphasize the qualitative aspects which are equal for all these systems. 



III. VIBRATIONS AND ATOMIC JUMPS IN THE GLASS 

The vibrations in a glass at low temperatures can be most easily studied in the harmonic 
approximation. As example we quench samples of up to 5488 atoms of a SSG to T = K. 
The atomic configuration will then correspond to a minimum of the potential energy. We can 
expand the energy in terms of the displacements from this minimum energy. The quadratic 
terms define a dynamical matrix 

D ij _ 1 d 2 E pot ({R}) 

where i,j denote the atoms and a, (3 the Cartesian coordinates. Diagonalization gives, in 
harmonic approximation, the frequencies of the eigenmodes of vibration and their eigenvec- 
tors, i.e. their spatial structure. Fig. 1 shows the vibrational spectrum, Z(u), together with 
the corresponding Debye spectrum calculated from the elastic constants. The area between 
the two curves is due to the excess low frequency vibrations, typical for glassy structures. 
In a plot of Z(y) / v 2 one finds a maximum around v = 0.1, the boson peak. The eigenmodes 
at the boson peak frequency have a complicated structure. At the lowest frequencies they 
can easily be identified as resonant (quasi-localized) vibrations. They can be decomposed 
into sound waves and soft local vibrations WQ. The latter are in full accord with the pre- 
dictions of the soft potential model ||. The interaction with the local vibrations causes an 
attenuation of the sound waves. With increasing frequency the Ioffe-Regel limit is reached, 
the phonon mean free path drops to the wavelength. The increasing interaction between the 
modes leads to a level repulsion and to Ziv) oc v. In the Z{y)jv 2 representation, or in the 
inelastic scattering intensity, this gives a drop oc v which causes the boson peak maximum. 

It is important that the cores of the resonant vibrations are local but extend over many 
atoms. Any single atom is stable against isolated displacements. The atoms in the cores of 
resonant modes participate also in high frequency localized vibrations. This is typical for 
resonant modes caused by stresses due to configurational "defects". They are not due to 
some atoms weakly coupled to the rest. 

The soft modes in glasses originate from some atomic configurations where some direc- 
tion in the 3A-dimensional configuration space is locally soft. This softness should not be 
confused with the one of the sound-waves which are soft because all atoms move in phase. 
In the extreme case, the group of atoms forming the center of the soft quasi-local vibration 
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is stabilized by the embedding matrix of the rest of the glass in a position lying between 
minima of the potential energy given by its near neighbors. To illustrate this, the dashed 
line in Fig. 2 shows the average potential energy of the 61 atoms which are most active in the 
given mode, {U(x)) coie . Atoms are considered as active in a given mode if their amplitude, 
| e n | , is at least 30% of the maximal atomic amplitude in the mode. The partial potential 
energy of these active atoms is indeed double- well shaped with minima at x m ~ ±1.3, which 
corresponds to maximal displacements of individual atoms by |R™ — Rq | « 0.2 — 0.3cx from 
the equilibrium configuration. This maximal atomic displacement is of the order of the one 
observed in local low temperature relaxations [IB|. 



At finite temperatures one observes these as aperiodic transitions from one local config- 
urational minimum into another. To visualize these we monitor the total displacement AR 
from a starting configuration and define 



Ai2(t)= /£(R"(t)-R"(0)) 2 (4) 

y n 

where R n (t) is the position vector of particle n at time t and R n (0) is the one at the 
starting or reference configuration. AR(t) oscillates due to the vibrations and changes due 
to relaxations, i.e. due to the transitions from one local energy minimum to another. An 
example is shown in Fig. 3 for the SSG at two temperatures. At T = 0.02T g the glass clearly 
jumps between three configurations. Let us denote the configurations A, B, C. The jump 
sequence is 5— >C^B^C^>A— > C — > A. By quenching to T = K we find the 
potential energy differences AE B - AE A = 1.74 x 10~ 3 e and AE C - AE A = 1.81 x 10" 3 e. 
These energy differences are of the order of the temperature. The corresponding spatial 
distances between the configurations, Eq. (H), are AR A - B = 1.63a, AR A ~ C = 0.96a, and 
AR B ~ C = 1.0a, i.e. they are of the order of the nearest neighbor distance Rnn — l-la. The 
maximal distance an individual atom travels in these jumps is only 0.3cr, about a quarter of 
Rnn ' Such relaxations can be observed experimentally e.g. as telegraph noise in the electric 
resistivity of point contacts [T7| . Increasing the temperature by a factor of four the average 



displacement in the initial configuration doubles as expected for a vibrational mean square 
displacement. The jumps seen at the lower temperature can no longer be resolved and new 
jumps over larger distances are observed. 

All observed relaxations are collective jumps localized to 10 or more atoms forming 



twisted chain-like structures with some side branching |16j . An example of such a structure 



is shown in Fig. 4 |L8| . The chain structure is a consequence of the dominance of the nearest 
neighbour bond in close packed structures. A low temperature relaxation is only possible if 
these bonds are not strongly compressed. Using the squared atomic displacement as "mass" 
we find a gyration radius of about 5 nearest neighbour distances for the chain of jumping 
atoms. The effective dimension is about 2 as one would expect for a twisted chain with 



side-branching. We observe similar structures in amorphous Se [B?jr7 and CuZr |T2 

In materials with a different inherent structure these relaxations will be different. In 
SiC>2 we would expect them to be collective twists of tetrahedra as has been postulated for 



the soft vibrations [21|. The origin will be, however, the same, namely local stresses leading 
to a softness in one direction of the multidimensional configuration space. 

Increasing the temperature we observe a marked increase of the number of atoms par- 



ticipating in a single jump 22 . 
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The local relaxations show two important correlations. First, in accord with the soft 
potential model they are strongly correlated with the soft vibrations [HJ16]. Secondly, there 



is also a strong but not full correlation between subsequent jumps in the same part of the 
sample. This latter correlation indicates that one jump can trigger another one which leads 
to a correlation in time, i. e. bursts of jumps ||23|| . It also leads to a slow increase in 



time of the number of atoms which have moved significantly, the dynamic heterogeneity. 
We illustrate it in Fig. 5 where we show the mobile atoms at T ~ 0.15T 9 after some ns 
(in Ar units). The effect is again independent of the specific material studied, whereas 
the size of the region at a given time and temperature will depend on the particulars of 
the inter-atomic interaction. The atoms which have moved in a given time interval form a 
complicated structure. From the time dependence of a suitably defined correlation function 
a fractal dimension has been derived for a binary soft sphere glass ||24j| . The time dependence 
of the dynamic heterogeneity is reflected also in the time dependence of the non-Gaussianity 
discussed further down. 



IV. SIMULATIONS IN THE LIQUID 

Computer simulations in the liquid state, especially with respect to MCT, have been 
discussed in numerous papers and reviews, see e.g. |TT|j25|j26[| . The predictions of MCT are 
at least semi-quantitatively reproduced in these calculations. We have not done extensive 
tests. However, quantities such as the intermediate self-scattering function and the self 
diffusion constant follow the trends of MCT for both our super-cooled CuZr [f27| and Se pg 
Tests of MCT are widely discussed elsewhere in this conference. 



In the liquid state it is no longer straightforward to separate single jumps. To obtain 
information on the atomic structure of the motion one can study the difference between 
structures some ps apart. To remove the effects of vibrations these configurations have to 
be averaged over a typical vibrational period. Fig. 6 shows an example for under-cooled 
liquid Cu33Zr 67 . There are two striking effects. First, clearly the smaller Cu atoms are much 
more mobile. Secondly we observe string like structures similar to the ones shown for the 
amorphous materials, Figs. 4 and 5. These structures indicate a high degree of collectivity 
in the motion also in the under-cooled liquid. They have been studied extensively for a 
binary Lennard- Jones system by Donati et al. . These authors report a marked increase 
of collectivity when the liquid is quenched to T c . These mobile strings again will lead to a 
dynamic heterogeneity. If one marks the atoms which have moved in a given time intervall 
by further than some cutoff distance, one observes structures similar to the one shown in 



Fig. 5 for the glassy state |30[. Experimental evidence for collective particle motion was 
gained by neutron scattering |T| its conclusiveness was, however, disputed [[32(1 . Additional 
evidence is again the vanishing isotope effect in diffusion 



V. ISOTOPE EFFECT IN DIFFUSION 

The isotope effect is the most direct experimental probe of collectivity. In a monatomic 
liquid the diffusion constant can be written 



D = D f(T,p) = D*f(T,p)/V^ (5) 
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where T is the temperature, p the atomic density and m is the mass of the diffusing particle. 
In the case of different components and different isotopes, considered here, the situation 
is more complicated. At low densities and high temperatures when diffusion is dominated 
by binary collisions the kinetic approximation should hold and Eq. [5] should apply approx- 
imately for each constituent. Lowering the temperature or increasing the density effects of 
collective motion will gain importance. 

A frequently used measure of this effect is the isotope effect parameter E K) 



Kb = - (6) 



where the index I denotes the different chemical components and a and (3 denote different 
isotopes. For a tracer atom with a mass of rrf + 5m e (average mass m e ) the change of the 
diffusion constant is then in linear approximation |jf 



An isotope effect parameter of rs 1 indicates uncorrelated single particle motion whereas 
collectivity results in E^ B — > 0. 

Formally one can describe the isotope effect by an effective mass 

K) eff = m e a + (N D - l)m (8) 

where No stands for the effective number of particles moving cooperatively. Such a definition 
can be justified well in the solid state. There, one has well defined diffusional jumps and 
No is defined by the multidimensional jump vector. Here we introduce No only formally. 
The mass dependence of the diffusion constant is then given by 



Di = f(T,p)/J(mi) eS . (9) 

From Eq. one gets then E l a3 » 1/N D . 

Evidence from earlier measurements of the isotope effect pl| is conflicting. Progress 
was made by simultaneously measuring the diffusion of the tracer atoms 57 Co and 60 Co |[37|| . 
Using this technique for diffusion of Co in amorphous C076.7Fe2Nb14.3B7 a value E = 0.1 was 
found indicating a high degree of collectivity. In contrast for self diffusion in crystalline Co 
one finds E = 0.7 . There, diffusion is by a vacancy mechanism which involves essentially 
single particle jumps with not too large displacements of the neighbors. The technique was 
also applied to a super-cooled melt of Zr4 6 .7Ti 8 . 3 Cu7.5NiioBe27.5, and again a very low isotope 
effect was observed |33 . 



We have calculated the isotope effect for a monatomic and binary Lennard- Jones 
liquids. Fig. 7 shows the isotope effect parameters for both components of a binary L J-liquid, 
consisting of 80% large and 20% small atoms as used in many other simulations, see e.g. 



1 1| • We find that the values are relatively low in the whole temperature range investigated, 
except for the high temperature ones of the smaller component. The behaviour for the larger 
atoms closely resembles the one observed in the mono-atomic LJ-system. These atoms form 
the backbone of the structure. The values for the smaller atoms show the same decrease 
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with temperature but are always clearly higher than those of the larger atoms. This reflects 
their smaller size. At the highest temperatures the small atoms seem to move through the 
matrix of larger atoms by binary collisions only. The observed high collectivity near T c is in 
agreement with the experimental findings of Ehmler et al. [33] who found E Pd 0.09 for Co 
diffusion in super-cooled liquid Zr46.7Ti8.3Cu7.5Ni10Be27.5- It is similar to the one found in the 



monatomic soft sphere glass at low temperatures [16|. The Decrease of the isotope effect with 



lower temperatures is proportional to the increase in density. The proportionality factors 
are, however, different for the two components. The effect of the density on the diffusion 
constant and on E are different, the dynamics clearly depends on density and temperature. 

The high collectivity both above and below the glass transition is related to the dynamic 
heterogeneity which we expect, therefore, to behave similarly. An experimentally accessible 
quantity which probes heterogeneity is the non-Gaussianity which we will discuss in the last 
section. 



VI. DYNAMIC HETEROGENEITY 

The change of dynamics upon cooling the melt towards the glass transition can clearly 
be seen in the self part of the van Hove function G s (r, t) which is related to the probability 
that an atom has moved by a distance r during a time t: 

P{r,t) = A7rr 2 G s {r,t) = (5 (r - |R n (t) - R n (0)|)) . (10) 

At high temperatures P(r,t) is, apart from the geometrical factor Anr 2 nearly perfectly 
Gaussian and broadens oc \ft. Upon cooling towards T c , and beyond, a tail to larger r-values 
grows with time. The self correlation function becomes markedly non-Gaussian and a tail 
to larger distances grows with time. Finally approaching T c additional structure evolves, 
particularly for the more mobile components, j|(],f|l],|ll],f4^,[27[] . As example Fig. 8 shows the 
distribution of the atomic displacements of Zr67Cu33 at T = 1000 K after a time t = 210 ps, 
i.e. during the early time of the so called a-relaxation. The second peak grows at a fixed, 
time independent, position roughly equal to the mean nearest neighbor distance. From 
this structure of P(r,t) one can conclude that there are preferred positions on the relevant 
time scales. This does, however, not necessarily mean that these positions are reached in 
a single jump, as frequently assumed. From an analysis of the evolution with time it has 
been concluded that the evolution of this structure in P(r,t) is accompanied with a strong 
increase of back correlation j|3) . 

This strong deviation of P(r, t) indicates that some atoms have been much more mobile 
than the average. This effect is called dynamic heterogeneity. To quantify it the non- 



Gaussianity parameter [44] is often used 



3<Ar 4 (t)> 1 

where < ... > denotes time averaging, Ar 2 (t) is the mean square displacement and Ar 4 (t) is 
the mean quartic displacement. This parameter is defined so that it is equal to zero when the 
atomic motion is homogeneous. Experimentally it can be obtained from the g-dependence 



of the Debye- Waller factor | 45| . It has also been calculated in numerous molecular dynamics 
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simulations of liquids, e.g. by Kob et al. for the binary Lennard- Jones system [|TTl , f46 
The limiting values for both times t = and t = oo is a 2 = 0. The latter limit reflects 
the ergodicity of the system for long times. Starting from t = 0, a^t) rises in general 
monotonically to a maximum from where it drops again monotonically. The maximum 
value is around 0.2 in the hot liquid and rises strongly in the under-cooled liquid, where a 
maximum value of 3 has been reported [[II] . The position of the maximum in time is at high 
temperatures in the ps range and correlates in the super-cooled liquid with the onset of the 
a-relaxation which is attributed to long range motion. This general behavior is observed in 
all systems we have studied. 

In the previous sections we have seen that on an atomic level the dynamics in the glass 
and the under-cooled liquids are very similar. To see this quantitatively for the heterogeneity 
we have calculated ^(t) for different temperatures above and below the glass transition. In 
Fig. 8 we show this for both components of the binary LJ-system in a log-log representation 
of t ■ a 2 (t) versus time. First we see the general trends discussed above. These hold not 
only in the liquid but also in the glass. Secondly we observe an asymptotic linear increase 
of ot2(t) oc y/i. For the lower temperatures this behaviour stretches over several ns (when 
one converts the L J- units to ones appropriate for Ar). This dependence can be explained by 
the above mentioned strong correlation between successive collective jumps |J47[| . One can 
envisage a chain of atoms jumping. This jump can either be reversed in a successive jump 
or might trigger a jump of another chain which will involve many but not all atoms which 
have jumped previously. 



VII. CONCLUSION 



We have shown by computer simulations of different materials that motion both in the 
glass and in the under-cooled liquid is highly collective. The dynamics in the glass at low 
temperatures can be described by the soft potential model which postulates similar struc- 
tures for tunneling centers, quasi localized vibrations and local relaxations. These local 
motions involve groups of ten and more atoms forming in dense packed systems predomi- 
nantly chain structures. These structures are closely related to resonant vibrations which 
cause the boson peak in the inelastic scattering intensity. Successive jumps in the same part 
of the sample are correlated. With increasing temperature the number of atoms jumping 
collectively increases in the amorphous state. Similar collectivity is also observed in the un- 
dercooked liquid where it increases upon quenching. The decrease of the isotope effect of the 
diffusion upon quenching to the glass transition indicates again an increase in collectivity, 
in agreement with experiment and other simulations. The collectivity of motion is related 
to the dynamic heterogeneity. The non-Gaussianity parameter, a measure of the dynamic 
heterogeneity, increases rapidly in the undercooled liquid and varies smoothly through the 
glass transition. 
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FIG. 1. Configurationally averaged vibrational density of states of the soft sphere glass (solid 
line) with n = 6 and Debye spectrum (dashed line) ||l4| 
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FIG. 2. Average potential energy of atoms in a single soft vibrational mode with frequency 
v = 0.0985, participation ratio p = 0.23. Solid line: potential energy averaged over all atoms in 
the system. Dashed line: partial potential energy, averaged over 61 atoms of the core of the mode. 



Dotted line: least squares fit of the partial potential energy by a soft potential polynomial | 15| . 
The average nearest neighbour distance is 1.1a. 
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FIG. 3. Total displacement relative to a local minimum configuration as function of time for 
one sample of the glass of Fig. 1 with N = 5488 at two temperatures. Please note the different 
scales for AR [16]. 




FIG. 4. Collective jump in the glass of Fig. 1 at T = 0.15T g . The initial positions of the atoms 
are shown by the dark spheres and the final ones by the shaded spheres. The bonds connect nearest 
neighbours. Shown are all atoms with more than 40% of the maximal atomic displacement fli~8f . 
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FIG. b. Initial positions in a sample of the glasses of atoms displaced by more than 0.3 of the 
average nearest neighbour distance during a time interval At = 3600(mo" 2 /e) 1//2 at T 
The total displacement is approximately 8 nearest neighbour distances [|l 



0.15T g . 




FIG. 6. Atoms in a melt of Cu33Zr67 at T = 1200 K with the largest displacements between 
time averaged configurations separated by 6.5 ps. Shown are Cu (light spheres) and Zr (dark 
spheres) atoms displaced by more than 1.6A and 1.45A, respectively [J22|. 
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FIG. 7. Isotope effect in a binary LS-liquid as function of temperature at equilibrium density. 
The diamonds and circles refer to the larger and smaller atoms, respectively. The arrow indicates 
the glass transition temperature p9f . 
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Distribution of the atomic displacements P(r, t) in Cu33Zrg7 near T c after 210 ps, 
The line shows a Gauss-fit to the small distance part 
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FIG. 9. Non-Gaussianity parameter multiplied by time for the binary LJ-system of Fig. 7 above 
and below the glass transition temperature. The solid and dashed lines refer to the larger majority 
and smaller minority components, respectively. The curves refer to the temperatures (from bottom 
to top): 0.88, 0.56, 0.48, 0.40, 0.36 and 0.32 in units of e/k. 
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